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ABSTRACT 

We report a detection of far ultraviolet absorption from the supernova remnant SNR 0057 - 7226 in the Small 
Magellanic Cloud (SMC). The absorption is seen in the Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum 
of the LBV/WR star HD 5980. Absorption from O VI A1032 and C III A977 is seen at a velocity of +300 km s'^ 
with respect to the Galactic absorption lines, +170 kms^^ with respect to the SMC absorption. The O VI A1038 
line is contaminated by H2 absorption, but is present. These lines are not seen in the FU SE spectrum of Sk80, 
only ~ 1' (--17 pc) away from HD 5980. No blue-shifted O VI A1032 absorption from the SNR is seen in the 
FUSE spectrum. The O VI A 1032 line in the SNR is well described by a Gaussian with FWHM = 75 km s"^ We 
find log N{0 VI) = 14.33-14.43, which is roughly 50% of the rest of the O VI column in the SMC (excluding the 
SNR) and greater than the O VI column in the Milky Way halo along this sight line. The A^(C IV)/iV(0 VI) ratio 
for the SNR absorption is in the range of 0T2-0.I7, similar to the value seen in the Milky Way disk, and lower 
than the halo value, supporting models in which SNRs produce the highly ionized gas close to the plane of the 
Galaxy, while other mechanisms occur in the halo. The A^(C IV)/A^(0 VI) ratio is also lower than the SMC ratio 
along this sight line, suggesting that other mechanisms contribute to the creation of the global hot ionized medium 
in the SMC. The O VI, C IV, and Si IV apparent column density profiles suggest the presence of a multi-phase 
shell followed by a region of higher temperature gas. 

Subject headings: supernova remnants — ultraviolet: ISM — stars: individual (HD 5980) 



1. introduction 

Supernovae are thought to be one of the main sources of 
the hot coronal gas in the interstellar medium (ISM) (McKee 
& Ostriker 1977; Sembach et al. 1997). There have been few 
absorption-line observations of the hot gas in individual super- 
nova remnants (SNRs), however, so global descriptions of the 
ISM have been forced to rely on models of evolving SNRs to 
reproduce the characteristics observed in absorption line studies 
of the Galactic disk and halo (Slavin & Cox 1992, 1993; Shull 
& Slavin 1994; Shelton 1998). In particular, the O VI AA 1032, 
1038 lines are very important probes of collisionally ionized 
gas near 3 x 10^ K. While several SNRs have been observed in 
O VI emission (Blair et al. 2000; Sankrit et al. 2001), few have 
been studied in O VI absorption. Jenkins, Wallerstein, & Silk 
(1976) detected O VI absorption in the Vela supernova remnant 
with the Copernicus satellite, but since then the opportunities to 
observe the spectral region containing the O VI lines have been 
limited. With the launch of the Far Ultraviolet Spectroscopic 
Explorer (FUSE) (Moos et al. 2000), this spectral window has 
been reopened. 

HD 5980 (Sk78, AV 229, I = 302°.07, b = -44°.95) is a 
Luminous Blue Variable / Wolf-Rayet star on the edge of the 
H II region NGC 346 (N66) in the Small Magellanic Cloud 
(SMC). De Boer & Savage (1980) first noticed an absorption 
system at Vis/? = +300 kms^^ arising in a highly ionized 
cloud in International Ultraviolet Explorer spectra of this star. 
Fitzpatrick & Savage (1983) suggested that the cloud might be 
an SNR moving toward HD 5980 at -150 km s^^. This in- 
terpretation was confirmed by a radio detection of an SNR in 
N66, SNR 0057 -7226 (Ye, Turtle, & Kennicutt 1991) coincid- 
ing with the X-ray source IKT18 (Inoue, Koyama, & Tanaka 
1983; Wang & Wu 1992). High velocity Ha emission was also 
seen at this position by Chu & Kennicutt (1988). 



Koenigsberger et al. (2001) observed HD 5980 with the 
Space Telescope Imaging Spectrograph (STIS) on the Hub- 
ble Space Telescope (HST). They detected the SNR in ultra- 
violet absorption lines, and were able to separate the absorp- 
tion into two components at V0 = +312 and +343 kms~^ 
(Visi? = +300 and +331 kms^^) in some of the lines. They 
also detected excess absorption at V© = +33 and +64 kms^^ 
(Visfl = +21 and +52 kms^^), suggesting that the high ve- 
locity gas is located on the far side of an expanding structure. 
The STIS spectrum contains absorption features of N V and 
C IV, which probe gas at —2 x 10'' K and ~1 x 10^ K, respec- 
tively. In this paper we present FUSE data on even hotter gas 
in SNR 0057 -7226. 

2. observations and data reduction 

HD 5980 was observed by FUSE on 2000 July 2. The data 
are archived in the Multi-Mission Archive at Space Telescope 
Science Institute as datasets P1030101 - P1030104. Spectra 
were taken through the large (LWRS; 30" x 30") apertures. 
Four individual exposures were combined for a total integra- 
tion time of 5734 s. A spectrum of Sk80 (AV 232), an 07 Iaf+ 
star (Walborn 1977) which lies close to HD 5980 on the sky, 
and which we use as a comparison star (dataset P10302), was 
observed on the same day as HD 5980 and reduced in the same 
manner. 

The FU SE instrument consists of four channels, two opti- 
mized for short ultraviolet (UV) wavelengths (SiCl and SiC2: 
905-1 100 A) and two optimized for longer UV wavelengths 
(LiFl and LiF2: 1000-1 187 A). The O VI Hnes are covered 
by all four channels, but the effective area is largest for data 
recorded in LiFl. It is the primary source of data considered 
here, although all the results were verified in LiF2 data. The 
C III line is covered by the SiC channels. We require that any 
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absorption feature be present in two channels to be considered 
real. 

The raw data were reduced using the standard FU SE caU- 
bration pipehne (CALFUSE vl .8.7) available at the Johns Hop- 
kins University. The velocity zero-point was set by shifting the 
Milky Way component of the H2 and low ionization hnes to 
kms~^, based on the results of Mallouris et al. (2001). The 
spectra were binned by three pixels (~6 kms~^ near O VI), 
and the nominal resolution of the data is ~20 km s~^ (FWHM). 
The relative wavelength solution is accurate to ^6 kms^^ on 
average. The signal to noise ratio in the binned data is ~ 50 
near the O VI lines. The continuum levels were chosen by fit- 
ting low order (<6) Legendre polynomials to Une-free regions 
of the spectrum near the lines of interest. 

3. RESULTS 

Figure 1 shows the observed absorption profiles of O VI and 
C m in the spectra of HD 5980 and Sk80. The two stars are 
separated by ^ 1', or ^17 pc, assuming a distance of 59 kpc to 
the SMC (Mathewson, Ford, & Visvanathan 1986). Absorption 
features of O VI A1032 and C m A977 are detected at +300 
km s~^ in the FUSE spectrum of HD 5980, but no similar fea- 
tures are seen in the Sk80 spectrum. The edge of the radio SNR 
extends very close to the projected position of Sk80 (Ye, Turtle, 
& Kennicutt 1991). Either the +300 kms~^ absorbing region 
does not extend as far as Sk80, or Sk80 is in front of the ab- 
sorbing gas. Absorption from hot gas in the Milky Way and the 
SMC is also seen in the spectra of both stars, and will be the 
subject of a future paper. TheO VI A1038 line at +300 kms"^ 
is contaminated by two H2 hnes, P(l)5-0 and R(2)5-0, but it is 
clearly much stronger in the HD 5980 spectrum than in Sk80. 

The Milky Way and SMC O VI A 1032 profiles in the 
HD 5980 and Sk80 spectra look very similar. Figure 2 shows 
an overlay of the two spectra, as well as the HD 5980 spectrum 
divided by the spectrum of Sk80 (the stellar continua near the 
O VI A 103 8 and C m A977 did not match as well). The Milky 
Way and SMC components divide out almost completely. The 
+300 kms~^ absorption toward HD 5980 remains, and is well 
described by a Gaussian with a FWHM of 75 kms~^. This 
corresponds to a temperature of ^2 x 10^ K if the line is broad- 
ened solely by thermal motions. We do not see any absorption 
from the approaching side of the SNR in the divided spectrum, 
identified in the STIS spectrum by Koenigsberger et al. (2001). 

Koenigsberger et al. (2001) detected two red-shifted com- 
ponents associated with the SNR in low ions, a strong com- 
ponent at Vlsr = +300 kms^^ and a weaker component 
at Ylsr = +331 kms-^ The absorption in the FUSE 
spectrum is located at Vlsr ~ +300 kms~^, correspond- 
ing to the strong component. We do not detect the Vlsr = 
+331 km subcomponent in C III or O VI. The absorption lines 
from this component seen by Koenigsberger et al. (2001) were 
very weak, so this component may be blended with the main ab- 
sorption system at +300 km s~ ^ in the FU SE data. Fitzpatrick 
& Savage (1983) detected a weak component at +325 kms"^, 
which if present in the FUSE data would be blended with the 
broad +300 kms"^ component. 

Table 1 Usts the measured equivalent widths and column den- 
sities for the O VI and C III lines, and upper limits on S VI and 
Fe II. No O I or H2 absorption is seen at +300 km s^^. The S III 
and S IV lines in the FUSE bandpass are blended with other 
absorption hnes, so we cannot determine whether any high ve- 
locity absorption is present. Colunon densities were calculated 



by integrating the apparent column density per unit velocity 
over the velocity range given, an approach which is valid if 
the intrinsic fine width is comparable to or broader than the 
instrumental line spread function (LSF) (Savage & Sembach 
1991). At temperatures of 3 x 10"'' K, where the abundance of 
O VI peaks in colhsional ionization equilibrium (Sutherland 
& Dopita 1993), the line width from thermal broadening (w 
30km s-i) is larger than the FUSE LSF width (« 20km s'^), 
so the condition is most likely satisfied. The C III abundance 
peaks at a lower temperature (T ~ 7 x 10^ K), and it is possi- 
ble that unresolved saturated absorption exists. In this case the 
measured column density of C III is a lower Umit to the true 
value. The uncertainty estimates in Table 1 include statistical 
noise fluctuations and modest continuum placement uncertainty 
(Sembach & Savage 1992). For the very faint N V line, con- 
tinuum placement may dominate the uncertainty, so its effects 
were evaluated by varying the continuum fit. 

We have reanalyzed the STIS data taken with the E140M 
grating from Koenigsberger et al. (2001) so that the equivalent 
widths could be measured over the same velocity ranges used 
for the FUSE data. These values are also given in Table 1 
(lines with A > 1200A were measured from STIS data). We 
also measured the O VI equivalent width and column density 
toward HD 5890 using the spectrum that had been divided by 
the Sk80 spectrum. Using the divided spectrum, we place an 
upper limit (3cr) of log A''(0 VI)< 13.43 on the approaching 
side, ~ 11% of the colunm density of the receding side. 

4. DISCUSSION 

Models of evolving SNR (Slavin & Cox 1992, 1993; Shel- 
ton 1998) predict the production of O VI in gas heated to 
T > 10'"^ K by shocks. The existence of O VI in the SNR 
requires shock speeds of at least ~ 160 — 170 kms"^ (Har- 
tigan, Raymond, & Hartmann 1987). If the absorption seen 
in the FUSE spectrum traces only the receding side of the 
shell, and the systemic velocity of the SNR is +176 kms"^ 
(Koenigsberger et al. 2001), then the expansion velocity is at 
least 124 kms"^, and may be higher if the sight fine is off- 
set from the center of the remnant. Chu & Kennicutt (1988) 
detected Ha emission separated from the SMC emission by 
+170 km s^b, so shock velocities this high certainly seem plau- 
sible. Assuming that A^(0 VI)/7V(0) < 0.2 (Sutherland & Do- 
pita 1993), and using log [O/H] = -3.85 for NGC 346 (Peim- 
bert, Peimbert, & Ruiz 2000), we find fliat iV(H+) > 8.1 x 10^^ 
cm~^. The 3-ct upper limit on O I from the STIS spec- 
trum is A^(0 I)< 3.0 X 10^^ cm^^, which corresponds to 
A''(H I)< 2.1 X 10" cm-2, so H+ZH" > 38. Koenigsberger 
et al. (2001) estimated NiH) = (4.3 - 12.0) x lO^^ cm'^ in 
the shell, close to our estimated value of iV(H+), so most of the 
gas in the shell may be ionized if the H" and H+ are cospatial. 

Figure 3 compares the O VI apparent column density profile 
with those of C IV, Si IV, Si m, and Si II convolved with a 
15 kms~b FWHM Gaussian to approximate the FUSE res- 
olution. The red sides of the profiles of the high ions have 
similar shapes, while on the blue side the higher ions extend 
over progressively larger velocity ranges, so that in terms of 
Yblue, the velocity of the blue edge of the profile, Vb;„e(0 VI) 
< Vb(„e(C IV) < YbiueiSi IV). The C IV has a distinct tail of 
absorption to negative relative velocities that extends almost as 
far as flie O VI absorption. Near +300 kms"\ the O VI/C IV 
ratio implies temperatures ~ (2.0 — 2.5) x 10^ K, while the 
C IV/Si IV ratio suggests T < 10^ K if the gas is in ionization 
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equilibrium (Sutherland & Dopita 1993). There appears to be 
gas at several different temperatures at this velocity, a conclu- 
sion supported by the presence of Si n, Si III, and Si IV. The 
C IV and O VI absorption at +250 km s~\ along with the ab- 
sence of absorption by lower ionization lines, suggests that this 
material is dominated by hot gas. A possible interpretation of 
this ionization structure is a multi-phase shell of swept-up ISM 
moving at the highest velocity (-1-300 km s~^), with the cavity 
behind it filled with hot gas. The C IV/Si IV ratio is ~5 in the 
shell, similar to that seen in the old SNR Radio Loop IV , and 
on the high side of the ISM average of 3.8 ± 1.9 (Sembach et 
al. 1997). The Une ratios imply temperatures lower than that 
indicated by the O VI line width, suggesting that non-thermal 
motions contribute to the line broadening. 

Table 2 lists the observed N(C IV)/N(0 VI) and 
NiC IV)/iV(N V) ratios for SNR 0057 -7226, for the SMC 
and Milky Way along this sight line, and the general disk and 
halo values. Also listed are the predictions of models of dif- 
ferent mechanisms for producing high ions. The observed 
NiC TV)/N(0 VI) ratio in SNR 0057 -7226 agrees well with 
the predicted ranges of models of evolving SNRs. The observed 
value of N(C IV)/7V(0 VI) for SNR 0057 -7226 is close to the 
Milky Way disk value, and very different from the halo value. 
The (C IV)/A(N V) ratio in the SNR also agrees with the disk 
value. While there may be a metaUicity effect on the observed 
ratios, the observed values for SNR 0057 - 7226 lend support 
to models of hot gas production (Shull & Slavin 1994) in which 
supemovae are largely responsible for the highly ionized gas 



in the Galactic disk, while different processes such as turbu- 
lent mixing layers (Slavin et al. 1993) and radiatively cooling 
fountain gas (Shapiro & Benjamin 1991) are responsible for 
the hot gas in the upper halo. The observed SNR ratios are also 
quite different from the general SMC ratios along this sight line, 
suggesting that mechanisms other than evolving SNR, such as 
those listed in Table 2, contribute to the extended hot gas in the 
SMC. 

The column density of O VI in the SNR is ~40-65% of 
the O VI column in the SMC component (excluding the SNR). 
Thus, only three such supernova remnants would be needed to 
expleiin all of the O VI in the SMC in this direction. Since 
this sight line is in the star-forming region NGC 346, it is quite 
plausible that multiple SNRs may exist along this sight line. 
The column density of O VI in the SNR is greater than that in 
the entire Galactic component in this direction. If SNR 0057- 
7226 is representative of supernova renmants in general, then 
they may be able to account for much of the O VI observed in 
the ISM of the Milky Way. 

This project benefitted from discussions with Charles Dan- 
forth. We thank Alex Fullerton and the FUSE Hot Star Team 
for reducing the FU SE data. This work is based on data ob- 
tained for the Guaranteed Time Team by the NASA-CNES- 
CSA FUSE mission operated by the Johns Hopkins University. 
Financial support has been provided by NASA contract NAS5- 
32985. 
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Fig. 1. — Observed line profiles for selected interstellar lines in the spectrum of HD 5980 and Sk80. The spectra have been normalized by dividing by the 
maximum flux value in the velocity range plotted. The velocities indicated by the vertical dashed lines are from left to right: The MUky Way at km s~^, the SMC 
at ~130kms-i, and SNR 0057 -7226 at ~300kms-^ AU three components of O VI A1037.62 absorption are blended with other lines, and the Milky Way and 
SMC components of C in A977.02 absorption are heavily saturated. 
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Fig. 2. — The O VI A1032 line in the spectra of HD 5980 and Sk80. The top panel shows the two continuum-normaUzed spectra, with HD 5980 as the soUd Une 
and Sk80 as the dashed Une. The bottom panel shows the HD 5980 spectrum after dividing it by the Sk80 spectrum. 
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Fig. 3.— A comparison of the O VI A1031.93, C IV A1548.20, Si IV A1393.76, Si III A1206.50, and Si H A1260.42 apparent column density profiles. The 
profiles have been scaled by the factors listed in the upper left corner. The STIS profiles (C IV, Si IV, Si HI, and Si II) have been convolved with a Gaussian with 
FWHM = 15 kms~^ to compare with the O VI profile observed by FUSE. The Si 11 and Si III lines are strong and may contain unresolved saturated structure, 
but are shown to demonstrate the velocity structure. The column densities for these lines should be considered lower limits. 



FUSE Observations of SNR 0057 - 7226 



7 



Table 1 



Line Strengths and Column Densities 



Line 


log /A- 


Velocity Range 




log A/'b 






(km s-i) 


(mA) 


(cm-2) 


O VIA103L93 


2.14 


230 - 350 


206 ± 12 


14.41 ±0.02 


O VI A1031 93'' 


2.14 


230 - 350 


191 ± 8 


14 36 ± 03 


VIA103L93^''* 


2.14 


0-150 


< 34 


< 13.43 


C III A977.02 


2.87 


260 - 360 


251 ± 5 


13.97 ±0.04 


S VI A933.38 


2.61 


260 - 360 


< 23 


< 12.83 


FeIlA1144.94 


2.08 


260 - 360 


< 14 


< 13.03 


FeniAl 122.52 


1.95 


260 - 360 


< 30 


< 13.53 


C IV A1548.20 


2.47 


230 - 350 


112 ±8 


13.53 ±0.03 


C IV A1550.77 


2.17 


230 - 350 


57 ±5 


13.49 ±0.04 


N V A 1242.80 


1.99 


230 - 350 


12 ±4 


13.03 ±0.12 


OIA1302.17 


1.80 


230 - 350 


< 24 


< 13.48 



'^The /-values are from Morton (1991) except for Fe n, which is from 
Howketal. (2000). 

''The uncertainties given are Icr estimates, and limits are 3cr estimates. 

''These quantities were measured on the HD 5980 spectrum after dividing 
by the Sk80 spectrum. 

''This is the Umit on absorption from the approaching side of the SNR. 



Table 2 
Column Density Ratios 



Location 


N{C iy)IN{0 VI) 


A''(C IV)/A^(N V) 


SNR 0057-7226 


0.12-0.17 


2.0^.5 


SMC (toward HD 5980) 


>0.71 


>5.1 


MW (toward HD 5980) 


-1.3^ 


~6.8'* 


MW (disk) 


0.08-0.26'' 


^3.8" 


MW (halo) 


0.46-1.91'^ 


4.3-7.0" 


Mixing Layer'^ 


1.0-3.4 


7.8-28 


Cooling Fountain'' 


0.1-0.5 


2.2-6.8 


Conductive Interface^ 


0.2-0.4 


1.8-2.6 


SNR models^ 


0.1-0.2 


1.9-2.3 



''C IV and N V from Sembach & Savage (1992). 

''From Spitzer(1996). 

"C IV and N V from Savage et al. (1997). 

''From Slavin et al. (1993), see Sembach et al. (1999) for details. 
"From Sembach et al. (1999). 

^From Borkowski et al. (1990), see Sembach et al. (1999) for details. 
sFrom Slavin & Cox (1992), see Sembach et al. (1997) for details. 



